The effect of the spin-orbit geometric phase on the spectrum of Aharonov-Bohm 
oscillations in a semiconductor mesoscopic ring 
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Taking into account the spin precession caused by the spin- 
orbit sphtting of the conduction band in semiconductor quan- 
tum wells, we have calculated the Fourier spectra of conduc- 
tance and state-density correlators in a 2D ring, in order to 
investigate the structure of the main peak corresponding to 
Aharonov-Bohm oscillations. In narrow rings the peak struc- 
ture is determined by the competition between the spin-orbit 
and the Zeeman couplings. The latter leads to a peak broad- 
ening, and produces the peak splitting in the state-density 
Fourier spectrum. We have found an oscillation of the peak 
intensity as a function of the spin-orbit coupling constant, and 
this effect of the quantum interference caused by the spin ge- 
ometric phase is destroyed with increasing Zeeman coupling. 
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The spin-orbit interaction (SOI) gives rise to a geomet- 
ric phase in the quantum amplitude of a particle prop- 
agating along a closed trajectory. In ideal ID rings this 
can lead to quantum oscillations of transport parameters 
similar to the Aharonov-Bohm effect g,0] . In disordered 
conductors the interference between two time reversed 
paths produces the oscillation of mean conductance anal- 
ogous to the Altshulcr-Aronov-Spivak effect. Such oscil- 
lation in 1-D systems was shown by Meir, Gefen and 
Entin-Wohlman p] , and in systems of higher dimensions 
by Mathur and Stone B. Besides, SOI also modifies 
the shape of the Aharonov-Bohm and Altshuler-Aronov- 
Spivak oscillations. For a disordered material, the mean 
conductance is obtained as an ensemble average over a 
large number of measurements on different samples. One 
can try to detect the quantum effects associated to the 
spin-orbit phase by measuring the oscillations of mean 
conductance when the spin-orbit coupling strength or the 
external magnetic field is varied. To our knowledge, such 
experiments have not yet established any evidence of the 
spin-orbit geometric phase. 

For a disordered material, if one takes the Fourier 
transform of the mean conductance {g{B)) as a function 
of the external magnetic field B, the spectrum is domi- 
nated by the Altshuler-Aronov-Spivak oscillations which 
is periodic in magnetic flux with a period hc/2e. On the 
other hand, if one takes first the Fourier transform 5(1^) 
of a measured conductance g{B), and then performs an 
ensemble average {\g{i^)\) of the Fourier amplitude, one 



would expect that the so-derived spectrum will exhibit 
a main peak corresponding to the Aharonov-Bohm oscil- 
lations with a period hc/e in magnetic flux B. Conse- 
quently, the average of Fourier amplitude, {\g{i')\) repre- 
sents correlations of conductances measured at different 
magnetic fields. The dependence of these correlations 
on the SOI can then manifest itself in the shape of the 
mean peak. In a recent experiment ||6| on mesoscopic 
rings made from a AlSb/InAs quantum well structure, 
the data were analyzed in this way for the first time, 
and a split of the main peak in the measured spectrum 
(Iff(^)l) was observed. The authors of Ref. g| have con- 
jectured that the observed splitting is due to the strong 
Rashba SOI in a doped AlSb/InAs quantum well. 

While the experimental results in Ref. [|6| remain to 
be explained, in our opinion, the fundamental question 
needed to be answered is how the Fourier spectrum of 
the conductance correlations, particularly the shape of 
the main Fourier peak, is influenced by the interplay 
between the spin-orbit phase and the external magnetic 
field. This is the aim of the present Letter. The more 
suitable starting point for such a theoretical analysis is 



(I5MI') 



dBdB'e^-'^^-^'\g{B)g{B')), (1) 



where the interval between —Bq and Bq covers the region 
in which g{B) is measured. This interval is much larger 
than both the period of the Aharonov-Bohm oscillations 
and the magnetic field correlation range |^ of mesoscopic 
fluctuations. In this Letter we will calculate (|l|) in the 
diffusion regime of a disordered 2D semiconductor ring 
of width w and radius R. Our theory explains the main 
physical mechanisms which determine the shape of the 
Fourier spectrum. For the state-density correlator which 
partly contributes to (1^) we predict a split of the main 
Fourier peak when effect of the Zeeman interaction is 
not completely suppressed by SOI and mesoscopic fluc- 
tuations. After we point out that this same feature ap- 
pears in an ideal ballistic ring, we conjecture later that 
our theory also provides the physical origin of main peak 
splitting observed in a chaotic ballistic ring H . 

Since it is the Rashba H term rather than the Dres- 
selhaus [M term which gives the major contribution to 
the SOI in InAs based QW uM, for simplicity in this 
Letter we will neglect the Dresselhaus term. We will 



use the standard perturbation theory which was apphed 
previously B to analyze the conductance correlations. 
In the framework of this theory the correlator in (|l]) is 
expressed via two-particle propagators (Cooperons and 
diffusons) where one of the particles propagates at the 
magnetic field B and the other at B'. If one neglects the 
Zeeman interaction with the external magnetic field, the 
Cooperon propagator is a function of B+B', while the 
diffuson depends only on AB=B-B' . Hence, in the vicin- 
ity of the main peak at i'—27t'^R^/^o, where ^o—hc/e is 
the flux quantum, the major contribution to M) is given 
by diffusons. If the Zeeman interaction is taken into ac- 
count, the Cooperon propagator is no longer a function 
of B+B' alone. Nevertheless, Cooperon's contribution 
is relatively small if the Zeeman energy q^.b^q/ttR^ is 
much less than the Thoulcss energy Et—DTi/R^. There- 
fore, in the integrand of (|l|) we will retain only the part 
of the correlator associated to diffusons. 

Each of the diffusons is a component of a matrix 
D(r, r ) with four spin indices representing the spin states 
of an electron and a hole. Following Ref. ^ one can show 
that the correlator in the integrand of (|l]) is proportional 
to 



/ (frSr'{Re (Tr[D(r, r')D(r', r)]) 
+2rr[D(r,r')D^(r,r')]}. 



(2) 



In the above equation the trace is taken separately over 
the electron and hole spin indices. It is convenient to 
express D in the representation of the total spin S of the 
electron- hole two-particle system [|l^. In this represen- 
tation D is a 4x4 matrix with components Dmn- The 
indices m and n can have the values -1,0,1 (for the z 
component of the triplet) and s (for the singlet). 

In A3B5 semiconductors the spin-orbit coupling has 
the form iJ^o^hk-s for an electron having spin s and 
quasi-momentum k. This SOI and the Zeeman interac- 
tion determine the spin dependence of D, which can be 
written as lllSl 



r((-iv • V + - V AA + hk • S)'),„gD(r, r') 
-hZD(r, r') + — D(r, r') ^ 6{r - r') , 



(3) 



where t^ is the elastic mean free scattering time, r the 
phase breaking time, and the notation (• • ■)ang is an an- 
gular average over the Fermi surface. The term ZD(r,r') 
in the above equation is due to the Zeeman interaction, 
and the nonzero components of Z are 
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(4) 



We see that if B^B' the matrix Z contains components 
which mix the singlet part and the triplet part of the 
diffuson. 



We choose the gauge such that for the field difference 
Ai? the vector potential is AA=Ai?rt/2, where t is a 
unit vector tangential to the ring. Since Hso—Q:{kxSy- 
kySx) when the Rashba term dominates, the boundary 
conditions at the inner and the outer radii of the ring are 



»^D(r,r') 



am 



(t.S)D(r,r') 



0. 



(5) 



In a narrow ring with the width w much less than the 
radius R, D(r, r ) varies slowly across the annulus. If the 
elastic mean free path I is much shorter than w, such slow 
variation can be treated perturbatively in the diffusion 
approximation [Q. Using the boundary conditions (|^), 
after averaging over r, (H) is reduced to an effectively ID 
equation, where D depends only on the azimuthal angles 
if, ip' . This function can be expressed in the form 



D(^,^') 



a^S^iv'- 



EMn(¥'') 
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(6) 



Making use of the rotation properties of the spin operator 

exp[iS,((^ - <^')](S • n)exp[tS,{^' - (p)] = (S ■ n') , 

where n^=cos ip' and n' —simp' , we arrive at a set of 4x4 
algebraic equations for the components of the matrices 
M„. From these equations one can derive the following 
subset of equations which contains only the -1,0,1 com- 
ponents of the matrices M„ 

S2)](A0)^M„ 
(7) 



A(/)-C(S-N)]'M„ + [a + p2Mo(l 
-M„ = 1 , 



-ipS,(0 + 0')M. 



t^Et 



where a^w^/AR^, MQ^[{n-/^(j)f+a{l^(l)f+l/T^ET]~^, 
C,'^=l+{Ram* /U)'^ , and p—gm*/kflra. N is a unit vector 
with 7V:,=<(l-C-2)i/2, Ny=n'y{l-C-^Y/^ and N, = l/C. 
The three dimensionless parameters a, C, and p deter- 
mine the shape of the main peak in {\g{y)\^) given by (|^). 
a describes the dephasing due to the penetration of the 
magnetic flux A0=Ai?7ri?^/<I>o into the annulus of the 
ring. We should remind the reader that to calculate a 
we have assumed diffusive propagation of particles in the 
radial direction. However, it is reasonable to believe that 
the dependence on flux of the form a(A(/))^, as appear- 
ing in (M), is also valid for rings with ballistic transport 
along radial direction. In this case a can be considered as 
a phenomenological parameter. The spin-orbit coupling 
constant C gives the spin-phase winding number after a 
particle has traversed a closed path along the ring. The 
parameter p is related to the Zeeman interaction. It de- 
termines the amount of mixing between the triplet and 
the singlet components of the diffusion propagator. This 
mixing appears in (0) in the form p^Mo(l-S^). Hence, 
although p is small, the effect of mixing is enhanced by 
the resonance of the singlet diffusion mode and the Sz=(i 
component of the triplet. However, with stronger spin- 
orbit coupling the system is driven out of the resonance 
due to the term C(S • N) in @j. 
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FIG. 1. Fourier spectra of state-density for various Zeeman 
coupling p and SOI strength C,: p=0.006, ^=1 (solid curve); 
p=0.0033, C=l (dot-dashed curve); p=0.006, C=2.5 (dahsed 
curve). The frequency is defined as u)=2-k-v^o/-kB? . 

After substituting (g) into (g) and carrying out the 
integration, we need to perform a numerical summation 
over n in order to obtain the conductance fluctuations 
and the state-density fluctuations which, according to 
Ref. 1 15 1, are given by the flrst term of (0). The sum- 
mation over n is from -Umax to rimaa;=50, which gives 
converging results. The magnetic field Bq is set at a 
value corresponding to a flux of 300$o- Within a reason- 
able range of material parameters our numerical results 
depend only weakly on the value of Bo. 
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FIG. 2. Fourier spectra of conductance: /9=0.006, ^=1 
(solid curve); p=0.0033, C, = l (dot-dashed curve); p=0.006, 
(^=2.5 (dahsed curve). The frequency is ijj=2-K-v'^Q/-KR^ . 

Figs. 1-3 were calculated with a=2-10^"' and I/t^pEt = 
10~^. For the convenience of presentation, the frequency 
in Figs. 1 and 2 is defined as Lu—2TT-i'^o/TrR^ which is di- 
mensionless. Fig. 1 shows the main Fourier peak of the 



state density correlator. The splitting of the peak (solid 
curve) for p=0.006 and (=1 is due to the resonance of 
diffusion modes. As the Zeeman coupling is reduced to 
p— 0.0033 (dot-dashed curve), or as the resonance is de- 
tuned by increasing the SOI to C=2.5 (dashed curve), the 
phenomenon of peak splitting disappears. We have also 
found vanishing of this splitting when the value of the 
parameter a is enhanced, corresponding to a decrease of 
the magnetic field correlation range of mesoscopic fluc- 
tuations. While the splitting is seen in the state-density 
peak, it is absent in the diffusion coefficient spectrum. In 
the parameter regime considered here, the contribution 
to the conductance correlator from the diffusion coeffi- 
cient spectrum is larger than that from the state-denstiy 
correlator. Consequently, there is no peak splitting in the 
Fourier spectrum of conductance oscillations as shown in 
Fig. 2. 
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FIG. 3. The intensity of the main Fourier peak as a func- 
tion of the SOI coupling: p=10~* (cicular dots); p=5-10^^ 
(triangles); p=10~^ (squares). 

Our theory predicts a quantum oscillation of the inten- 
sity of the main peak in the Fourier spectrum of conduc- 
tance oscillation, as a function of the SOI strength. To 
demonstrate that its origin lies in the geometric phase, 
let us first set the Zeeman interaction p=0. In this case 
(0) can be easily diagonalized by choosing the spin quan- 
tization axis along N. The eigenvalues of CS • N are (, 
and — C. Including the singlet state, the eigenvalue 
is doubly degenerate. The dependence of {\g{i^)\'^) in (|l|) 
from ( can then be readily calculated, because the ±C 
can be absorbed by a shift of A(j} in corresponding com- 
ponents of the correlator. If we ignore the small change 
produced by this shift in the a(A0)^ term, after taking 
trace, {\g{j^)\'^) is found to be proportional to l-|-cos[(27r- 
uj)C]- At the center of the peak, uj=Q and so {\g{i^)\'^) 
oscillates with ^ as cos^ tt^. Since the Zeeman interac- 
tion breaks the time inversion symmetry and thus leads 
to an additional dephasing, the magnetic field correlation 



length is reduced and hence the dependence on the geo- 
metric phase gets weaker. As a result, the oscillating am- 
plitude of {\g{i')\^) decreases rapidly with increasing Zee- 
man interaction. The numerical result of this oscillation 
is shown in Fig. 3 for p=10~^ (cicular dots), p—b-lQ^'^ 
(triangles), and p=lQ^'^ (squares). 

The splitting of the main Fourier peak has been ob- 
served in a recent experiment on AlSb/InAs/AlSb quan- 
tum well rings with radii about 1 /ini p]. The material 
parameters for these samples are Z=10~'*cm, g—U, and 
a~10-9 eV cm Q, which give p^lQ-^ and ( =^3. With 
these values of C and p, our calculation does not yield 
the splitting of the main Fourier peak in both the state- 
density and the conductance. However, this is not a to- 
tal surprise because the samples used in Ref. |g] are very 
small with Z~i?, which is nearly in ballistic regime instead 
of in diffusion regime. Furthermore, in these samples the 
spin precession length h /am* is shorter than the elastic 
mean free path, for which our perturbative treatment of 
SOI is not valid. 

In order to judge how relevant are our qualitative re- 
sults to the above mentioned experiment, let us consider 
the opposite limit of an ideal ID ring. The period of 
the Aharonov-Bohm (AB) oscillations changes due to the 
dependence of the geometric phase and dynamic spin- 
phase on the magnetic field. We will consider first the 
effect of geometric phase. In the region of our interest 
p'(t><^{C'^-iy^'^ with p' defined as p'=gm*/kfRm, for an 
ideal ring one gets from Refs. pU2] the geometric phase 
9g=±TTp'{(^-l)^^/^<j)+C , where C is a constant indepen- 
dent of the magnetic field, and the ± signs refer to the 
two electron spin orientations. Combining the geometric 
phase to the AB phase 2tt(J), we see that SOI leads to 
a split of the AB oscillations in the transmittance of the 
ring into two oscillations with close frequencies. However, 
this frequency splitting is about two orders of magnitude 
less than the observed value [^. Next, we consider the 
split of the AB oscillation frequency caused by dynamic 
spin phases, which depend on the magnetic field in the 
form ±TT6kfR, where iSkf are the shifts of the Fermi 
wave-vector for up- and down-spin electrons. The shift 
irSkfR has its maximum value irp' (j)/2 in the absence of 
the SOI, and decreases with increasing SOI strength [|l],||. 
Even at the largest value 7rp'0/2, the corresponding split 
in AB oscillation frequency is of the order p' which is too 
small to explain the experimental value. 

The general feature of the peak splitting for the ideal 
ID ballistic ring is then similar to that for the diffusive 
2D ring. Furthermore, in the 2D ring the amount of 
the AB oscillations splitting as seen in the state-density 
correlations also decreases with increasing SOI. If we use 
the same value for p' in a ID ballistic ring and for p in a 
2D diffusive ring, which means the same strength of the 
Zeeman coupling in both systems, in the absence of SOI, 
the peak splitting in state density correlations is larger 
than that in the transmittance of a ID ring. This is due 



to the longer paths traversed by a diffusing particle in 
its random walk along a ring, and so acquiring a larger 
dynamic phase. We have reached the conclusion that 
the two quite distinct limiting cases have led to the same 
qualitative picture. Consequently, we conjecture that the 
peak splitting observed in near ballistic 2D samples |q| 
is due to the Zeeman interaction which is not completely 
suppressed by the SOI and the mesoscopic fluctuations. 
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